NASA Technical Memorandum 101574 


AEROSERV OEL ASTIC MODELING AND 
APPLICATIONS USING MINIMUM-STATE 
APPROXIMATIONS OF THE UNSTEADY 
AERODYNAMICS 


Sherwood H. Tiffany and Mordechay Karpel 


APRIL 1989 


(NASA-TH-1Qt574) A EROS EBV OEL ASTIC HO DEL IMG N89-24308 

AMD APPLICATION 0SIHG MINIMUM- STATE 

APPfiOXI HATIONS OF THE UNSTEADY AERODYNAMICS 

(NASA. Langley Research Center) 13 p Unclas 

CSCL 01C G3/05 0212638 


iw\s/\ 

National Aeronautics and 
Space Administration 

Langley Research Center 

Hampton, Virginia 23665-5225 


AEROSERVOELASTIC MODELING AND APPLICATIONS USING MINIMUM-STATE 

APPROXIMATIONS OF 
THE UNSTEADY AERODYNAMICS 


Sherwood H. Tiffany t 
Aeroservoelasticity Branch 
Structural Dynamics Division 
NASA Langley Research Center 
Hampton, Virginia 23665-5225 

and 

Mordcchay Karpel 
Senior Research Associate 
Department of Aeronautical Engineering 
Technion - Israel Institute of Technology 
Haifa, Israel 


Abstract 

Various control analysis, design, and simulation 
techniques for aeroelastic applications require the equations 
of motion to be cast in a linear time-invariant state-space 
form. Unsteady aerodynamics forces have to be 
approximated as rational functions of the Laplace variable 
in order to put them in this framework. For the 
Minimum-State method, the number of augmenting states 
representing the unsteady aerodynamics is a function only 
of the number of denominator roots in the rational 
approximation. Results are shown of applying various 
approximation enhancements (including optimization, 
frequency dependent weighting of the tabular data, and 
constraint selection) with the Minimum-State formulation 
to the Active Flexible Wing wind-tunnel model. The 
results demonstrate that good models can be developed 
which have an order of magnitude fewer augmenting 
aerodynamic equations than more traditional approaches. 
This reduction facilitates the design of lower order control 
systems, analysis of control system performance, and near 
real-time simulation of aeroservoelastic phenomena. 

Introduction 

The equations of motion of a flexible aircraft 
contain unsteady generalized aerodynamic force terms 
which are transcendental functions. The availability of 
efficient linear systems algorithms used in aeroservoelastic 
analysis and design has provided strong motivation to 
approximate the unsteady aerodynamic forces as rational 
functions of the Laplace variable, refs. 1-5. Such rational 
function approximations (RFA's) allow the 
aeroservoelastic equations of motion to be cast in a linear 
time invariant (LTI) state-space form, albeit with increased 
size of the state vector due to the RFA's. This increase in 
the number of states due to the RFA's is referred to as the 
aerodynamic dimension. There is always a trade off 
between how well the RFA's approximate the 
aerodynamic forces and the desire to keep the aerodynamic 
dimension small. The RFA formulations in the literature 
(e.g., refs. 1-12) have varying capabilities to perform 


such a trade-off. 

Currently there are three basic formulations used 
in approximating unsteady generalized aerodynamic forces 
for arbitrary motion using rational functions: 

1 . Least-squares (LS) - references 1 and 5 

2. Modified matrix-Pade' (MMP) - reference 
2, 3, 6, and 7 

3. Minimum-State (MS) - 8 and 12 

Extensions to these approaches were developed (refs. 9 - 
11) which included the capability to enforce selected 
equality constraints on the RFA's and of optimizing the 
denominator coefficients in the rational functions using 
nonlinear programming techniques. For the Minimum- 
State method, the number of augmenting states required to 
represent the unsteady aerodynamics is a function only of 
the number of denominator roots in the rational 
approximation; there are no multiple aerodynamic roots 
as there are for the Least-Squares formulation. Methods to 
determine critical frequency ranges based upon physical 
properties and to determine weightings of the individual 
tabular errors to improve approximations in those critical 
regions have been explored in reference 12. 

The evolution of advanced fighters has required 
that the aerodynamics, control systems, and structures 
disciplines be integrated into a unified aeroservoelastic 
technology which must be evaluated by sophisticated 
analytical methods and validated through the testing of 
wind-tunnel models. The Active Flexible Wing (AFW) 
program was a joint Rockwell/Air Force/NASA program 
(ref. 13) that involved the design, construction, and testing 
of a scaled aeroelastic wind-tunnel model. A follow-on 
cooperative effort between NASA and Rockwell is 
currently underway to investigate and evaluate the 
analysis, design and test methodologies associated with 
mulli-input/multi-output digital control laws. The 
purpose of this paper is to demonstrate the various 
Minimum-State aerodynamic modeling techniques by 
applying them to a current AFW wind tunnel model to 
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develop valid low order mathematical models for the 
control system design, evaluation and simulation tasks 
associated with the cooperative effort. 

The following pages present a review of the 
techniques involved. Results are also shown which 
demonstrate that good models can be developed which 
introduce an order of magnitude fewer aerodynamic state 
equations than more traditional approaches. 


tion of the errors, £ij(ik n ) = I Qij(ik n ) - Qij(ik n ) I, between 
the approximating function at p = ik n and the tabular 
value Qij(ik n ). It is possible that certain equality 
constraints or weighted constraints might be desirable to 
impose on some of the tabular values (such as the steady 
state ( k = 0) point, shown in the shaded box). 

Review of Rational Function Approximations to 
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Unsteady Aerodynamic Approximations 


In order to express the equations of motion in 
first-order form, the unsteady generalized aerodynamic 
force (gaf) coefficients (Qij) must be defined in the 
complex Laplace ^-domain, or nondimensionalized Laplace 
p-domain, p-(c!2u)s\ here u is the freeslream velocity and 
c is a reference length (typically the mean aerodynamic 
chord). Since the gaf coefficients are computed at 
specified values of reduced frequencies k n „ it is necessary 
to generate complex ^-domain functions by approximating 
the tabular values, Qij(ik n ), of the gafs as closely as 
possible, and then employing the concept of analytic 
continuation in a region near the imaginary axis. Figure 
1 depicts the approximating process for a single element 
Qij of Q , where 


k 

Qij(ik n ) 

Qij(ik) 

Qij(ik n ) 


eij(ikn) 


reduced frequency, (cl2u)a> 
reduced-frequency domain tabular data 
(identified by the open circles) 
approximating curve, Qij(p) , for 
p = ik (corresponding to the solid line) 
points along the approximating curve at 
reduced frequencies, k n , corresponding to the 
tabular data (solid dots) 
is the approximation error between two 
corresponding points (denoted by an arrow 
between the points). 
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Figure 1 Approximating aerodynamic tabular data for a 
single element of Q, employing constraints 


a 

The approximating function Qy(p) is determined 
in such a way as to minimize some least-square combina- 


The most common form of the approximating 
functions used currently for each generalized force 
coefficient, Qy , is a rational function of the nondimen- 
sional Laplace variable p . The normal form chosen is 
one in which the numerator polynomial is order 2 more 
than the denominator. This gives rise to the following 
partial fraction form for each element, Qy , where ni is 
the number of partial fractions and corresponds to the order 
of the overall denominator polynomial: 


n L 

+ +2) ^ p + bl 


(la) 


which can be rewritten in the Laplace domain as 


n L 

+ j^S A U + 2)^j 


where 


s + 


b n = ir b £ 

\ =A o 

A = ~A, 

I 2 u l 

A 2 = (TuJ A 2 

A l + 2 = A £+2 


(lb) 


Because tabular data are determined for specified 
values of reduced frequency, k n „ the Qy are actually 
defined only for these values of the nondimensionalized 
Laplace variable p = ik n . It is desirable to use the above 
rational function form for the approximating functions in 
order to convert the transcendental equations of motion 
into linear, time-invariant, state-space form. The fact that 
the numerator is only order 2 more than the denominator 
implies the number of resulting equations associated with 
the aerodynamics is a function of the order of the 
denominator polynomials. The Ao, Ay, and A2 terms can 
be included with the stiffness, damping, and mass terms, 
without any additional equations. 
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^ When the matrix equations, which combine all 
the Qij 's, are formulated, the character of the partial 
fractions are expressed in various ways (ref. 11). The 
first, 


A - a 

Q (p) = A 0 + Aj p+ A 2 p 2 + Q l , where 


Q, - I A t 


iL £% ^ + 2 ) P + be 


( 2 ) 


is referred to simply as the "Least-Squares" (LS) 
formulation in which the number and values of 
denominator coefficients, bj>, are fixed for all Qij . The 
second. 


Qj(p ) 

where 


{A 0 } y + {A j p +{A 2 } y p 2 +Q Lj 

q _ y r A -1 p ( 3 ) 

Ql > ' ti ^ u + 2) tj p + bt i 


R = 


, £ =1,...,N 


(5) 


The resulting number of aerodynamic states 
introduced into the first order equations as a result of the 
MS formulation is N. It is usually a little larger than ni, 
but generally much smaller than either the aerodynamic 
dimension for the least-squares method, or that of 

the modified matrix-Pade', 


n E +n g 

I n L 
7=1 J 


The character of each Ql and the aerodynamic dimensions 
for the different formulations is summarized in table 1. 


for each column y, is a Modified Matrix-Pade 1 (MMP) 
formulation, which allows the number of denominator 
coefficients and their values to vary between columns. 
They are only fixed per column. The Minimum-State 
(MS) formulation, 

A o A 

Q (p) = A 0 + Aj p+ A 2 p 2 +Q l , where 

A 

Q l = D(pl -R) Ep (4) 

fixes the denominator coefficients as in the least squares, 
but additionally, the partial fraction numerator coefficients 
are determined as a coupled product of a premultiplying 
matrix, D, and a postmultiplying matrix, E. The diagonal 
matrix of aerodynamic roots is denoted by R\ i.e.. 


Constraints and Lag Coefficient Optimization 

Extended versions of the LS, MMP, and MS 
methods to fitting functions to aerodynamic tabular data 
are included in the ISAC program (ref. 14). The 
extensions include capabilities to select equality 
constraints and to perform a nonlinear optimization of the 
denominator coefficients (refs. 9-11). The equality 
constraints allow for more realistic modeling of the 
aerodynamics and for improved fits at critical points (such 
as flutter). Nonlinear optimization of the denominator 
coefficients allows improvement in the approximations 
without increasing their number by determining a better 
set of coefficients than might be chosen apriori. 


Aer odynamic Method 
Least-Squares 


Modified Matrix Pade' 


Minimum-State 


Table 1 RFA Matrix Formulations 
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Character of Ql 

Common denominator coefficients in each 
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Qij 


Different number of and values for denominator 
coefficients for each column , Qj 
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Common denominator coefficients in each Qij 



« D n E ljP 

t, 


Aero dynamic Dim e nsion 


n 


4 


• n 


L 


"4 


j 


+ no 

I n L 

= 1 J 


N 


3 


Minimum-State. -Equations of Motion 


The Minimum-State matrix equation for the 
unsteady aerodynamics is given by equation 4. The 
resulting open-loop state-space equations of motion are: 
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where x a =[i/ - ^t~r\ 

The 4 represents structural modes and the 8 represents 
control modes. The submatrices are defined as follows: 

* - « - <- t ) 2 [* 2 ] 

K -K 

‘ r -' Hl , 

0 ‘ G 


(mass matrix) 
(stiffness matrix) 

(damping matrix) 


<7 


'8 

nondimensional gust velocities 
dynamic pressure 


D , E , and R were defined in equations (4) and (5). The 
subscripts refer to column subsets involved; £ refers to 
vibration modes, 8 refers to control modes, and g to gust 
modes. 

Phy si cal Wei gh tin g 

In ISAC (ref. 14) the inclusion of constraints is 
optional and those selected are imposed on a per column 
basis. A recently developed program, MIST (Minimum 
State), which has been interfaced with the ISAC program, 
provides the capability to include weights, 


w .. - w .. 

ijn ijn 

where 
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*0 =ma n X {\ Q ij( ik n) f'ijn } and 0 <w cut < 1 
on each £[j(ik n ) to weight the importance of the fits at 


each tabular point Qij(ik n ). Each Wy n is a measure of 
importance which allows the fits to be improved at some 
points (at the possible expense of others) based upon 
physical properties without actually enforcing equality at 
the specified points (ref. 12) The method developed 
therein determines the measure of importance based upon 
partial derivatives of selected open-loop parameters with 
respect to the weighted term at a specified design point. 
For the vibration modes the weight of the fit at each value 
of k is determined by the effective influence on the system 
determinant; for control modes, by the effect on Nyquist 
gains; for gust modes, by the effect on the response to 
continuous gusts; and for hinge moment terms, by the 
hinge moment response to control surface or gust 
excitations. The physical weights are ^-dependent, but 
they are independent of normalization of the generalized 
aerodynamic force coefficients. The w cut in equation (7) 
provides a means to specify the minimum of the 
maximum weighted magnitudes 

of each Qij . It can be shown that 

w cut ^Q'ij . 

Hence, a value of w cu t = 0 allows the full influence of 
physical weighting by setting the minimum = 0; whereas 
w cut = 1 forces the most sensitive point for each Qij to 
be normalized or have a weighted magnitude of 1. 

The current MIST formulation requires the 
imposition of three equality constraints on each element 
Qij. These equality constraints can be selected from 
various options, but must be the same for all Qij . For 
the results shown herein, this requirement is not severe. 

Data normalization is the basic measure of error 
used in ref. 9-11 which enforces fits equally (except for 
equality constraints) over a selected frequency range of 
tabular values in order to best meet the criterion for 
analytic continuation. However, since the concept of 
analytic continuation used to extend frequency domain 
functions, descibed by a finite set of tabular aerodynamic 
data, is approximate at best, physical weighting with 
optimization provides another tool for extending these 
functions in some reasonable fashion. Furthermore, 
varying the minimum of the maximum weighted 
magnitude, by varying the value of wcut between 0 and 1, 
provides a mechanism which allows the approximations 
to incorporate some combination of physical weighting 
and simple data normalization. 

AFW State-Space Modeling 

The AFW project requires various first-order 
models for different purposes such as low-order models for 
control law design and near real-time simulation and 
higher-order models for control system evaluation in batch 
simulation. 
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For both types of simulation it was desired to 
generate a reliable mathematical model in which the mass 
matrix was constant, i.e. had no aerodynamic terms, in 
order to avoid repetitive matrix inversion. This implied 
that the A 2 coefficients in the rational approximations had 
to be zero for all rigid and elastic modes. To meet the low 
order requirements of design models and near real-time 
simulation, MS models were generated using 2 
denominator coefficients. This resulted in an 
antisymmetric aeroservoelastic model with a total of 38 
states. These states consisted of 22 rigid and elastic mode 
displacement and rate states, 12 actuator states, 2 gust 
states, and 2 aerodynamic states. 

Discussion of Results 

The following results are for the antisymmetric 
AFW model at Mach .9. There are 1 1 vibration modes 
(the roll mode is included in this group), 4 control modes 
and one gust. In each case, the aerodynamic 
approximations arc constrained to match at k= 0, the 
imaginary parts are are matched at k- 0.005, and the A 2 = 
0. Furthermore, there is no constraint on the minimum 
that the maximum weighted magnitudes of Qij must be 
(i.e., w cut = 0.0), corresponding to full inclusion of 
physical weighting. 

Figure 2 is a set of three dimensional plots (from 
two different views) of the physically weighted 
magnitudes of each Qij(ik n ) These weights were 
determined for a value of dynamic pressure (q ) near flutter. 
On each plot, Qij> (j=7 to 11) are shown in a group for 
each j. For the vibration modes, j increments from 1 to 
77, corresponding to the generalized coordinates. For the 
control modes,; increments from 72 to75, corresponding 
to the leading edge inboard (LEI) , trailing edge inboard 
(TEI), leading edge outboard (LEO), and the trailing edge 
outboard (TEO) control surfaces. There is one gust mode, 
corresponding to ;=76. The peaks correspond to the 
tabular values of reduced frequency of highest influence on 
system characteristics. For example, for the roll mode 
(column 1) the system determinant is most sensitive to 


changes in the generalized aerodynamic force coefficient 
near the third k value, k= 0.01, in the first row; i.e., the 
system is most sensitive to errors in Qj](k 3 ). 

Some of the rows have nearly zero weighted 
magnitudes for all values of k , which indicates that the 
system has minimal sensitivity to variation in these 
elements. The view on the left indicates that for the 
vibration modes, the diagonal terms, Qu , are the most 
sensitive. Referring to Table 2 for the values of in 
vacuum natural frequencies (ov)» the view on the right 
indicates that the area of sensitivity is at k values near k v 
of each mode. This can be seen more easily in figure 3, 
which shows the weighted magnitudes of the most 
sensitive vibrational mode elements, Qu , at each tabular 
value of k. According to this plot, the system is most 
sensitive to errors in Qss at reduced frequencies between 
0.6 and 0.8. This plot also indicates that the generalized 
forces probably should have been calculated for additional 
values of reduced frequency between 0.6 and 1.5, as the 
peaks are not clearly defined for some of the modes in this 
region. As indicated by figure 2, the same denominator 
coefficients based on sensitivity for modes 1 through 5 
should work for the control modes and gust in for this 
model. They have the same sensitivity range as the first 5 
modes; i.e., they are not very sensitive to fits past k= 0.6. 
In order to obtain a reliable model for all the modes 
selected however, it appears that the range of tabular 
values (0.0<£<2.0) was appropriate for this configuration. 

The determination of an optimal set of 
denominator coefficients for a Least-Squares fit is a fairly 
quick process (reference 11). For the Minimum-State 
approximation, however, nonlinear optimization of the 
denominator coefficients requires a three-fold, time- 
consuming iteration process. In lieu of using nonlinear 
optimization, figure 3 was used to help select denominator 
coefficients for several cases. The selections, although not 
optimal, produced good fits, and results are presented in 
table 3 and figure 4 in order to show the reader the type of 
results which can be obtained with various selections of 
denominator coefficients. 


Table 2 Frequency parameters (c o v ,k v ) at q = 0 and (cof,kf) at flutter 



(O v 

(rad/sec) 

k v 

*f' 

(rad/sec) 

*/ 

Most sensitive 
range of kt 

1 

2.39 

0.009 

9.34 

0.034 

[0.0,0.11 

2 

44.13 

0.159 

45.37 

0.164 

[0.1.0.2] 

3 

49.37 

0.178 

53.93 

0.195 

near 0.2 

*4 

82.23 

0.297 

73.05 

0.263 

[0.2,0.3] 

5 

101.53 

0.366 

101.64 

0.367 

near 0.4 

6 

173.04 

0.625 

200.11 

0.722 

[0.6,0.81 

7 

240.65 

0.869 

217.93 

0.787 

near 0.8 

8 

248.93 

0.899 

248.79 

0.898 

[0.8,1.0] 

9 

258.77 

0.934 

252.54 

0.911 

[0.9,1.01 

10 

313.74 

1.132 

313.95 

1.133 

[1.0,1.2] 

11 

326.88 

1.180 

348.58 

1.258 

near 1 .2 


* Flutter mode 

trange of k for each mode in which system determinant is most sensitive to error 
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Figure 3.- Weighted magnitudes of most sensitive aerodynamic elements 


Table 3 shows the choices made for the order of 
approximation equal to 1, 2, and 4 (equals the number of 
denominator coefficients). Also listed are the method for 
selection and the % errors in flutter q and frequency. TJjte 
word ’unconstrained' refers to k values for which the Q t j 
are not constrained; i.e., those other than £=0 and 
£=0.005. The words 'estimated most sensitive Qq ' refers 
to a case in which the estimated peaks (as identified in 
figure 3) are used to determine the most sensitive Qij . 
The k's corresponding to the Qij with 'estimated peaks' 
are also estimated. 


Figure 4 shows the open loop root loci of a 
baseline configuration and the best overall of the two and 
four-aerodynamic state cases as indicated in table 3. The 
baseline used the tabular, frequency domain generalized 
aerodynamic forces and was generated using a p-k 
determinant iteration process employing interpolated 
values (ref. 15). The other two were based on linear time 
invariant state-space methods. The Hassig form of 
implementation of the p-k method (ref. 16) becomes less 
accurate as the approximation Q(p) = Q(0+ik) degrades; 
hence the baseline result for the roll mode in figure 4 
should not be regarded as accurate. 


Table 3 

Flutter parameters for 

various minimum-state physically weighted approximations 

Order of 

Values of Selected 


<\f 

(Of 






1 

0.8 

(unconstrained) k of most sensitive Qij 

0.06 

0.27 


1.0 

Mid-range k 

-0.07 

0.29 


0.722*t 

kf of most sensitive Qij 

0.13 

0.26 


0.625 

k v of most sensitive Qij 

0.24 

0.25 

2 

0.2, 0.8* 

(unconstrained) k's of most sensitive Q { j s 

-0.07 

0.32 


0.625, 0.869 

k v 's of most sensitive Qij s 

-0.87 

0.38 


0.297, 0.625 

(unconstrained) k v 's of most sensitive Qij's, for different j 

0.32 

0.12 


0.625, 1.180 

k v ’s of estimated most sensitive Qij s 

-0.78 

0.37 


0.4, 1.5 

Mid-range k's of least sensitive Qij's 

-0.55 

0.36 

4 

0.2, 0.6, 0.8, 1.0 

(unconstrained) k's of most sensitive Qij's 

0.64 

0.00 


0.2, 0.8, 1.0, 1.5* 

(unconstrained) k's of most sensitive Qij's, for different) 

1.00 

0.05 


0.2, 0.73, 0.8, 1.22 

k's of estimated most sensitive (unconstrained) Qij's 

0.59 

0.01 


0.009, 0.625, 0.869, 1.180 

k v 's of estimated most sensitive Qij s 

0.09 

0.12 


0.159, 0.366, 0.899, 1.132 

k v 's of least sensitive Qa's 

0.98 

0.05 


*Best frequency and damping properties at flutter for all modes for corresponding order of approximation 
( Very little difference overall between this and other three cases 
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Imag 

(rad/sec) 



(0 < q < 400 psf) 

%difference in flutter q < 1 .0 (=3.25 psf) 
^difference in flutter frequency < .4 (=.3 rad/sec) 


Figure 4.- Stability analysis root loci for the AFW antisymmetric modes at Mach 0.9. 
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Table 4 Comparison of optimised physically weighted Minimum-State approximations 


Order 

Values of Selected 
Denominator Coefficients. 

Physically Weighted 

Method of Selection Total Error 

Normalized 
Total Error 

m 


1 

1.225 

Nonlinear Optimization 

0.32 

17.75 

-0.18 

0.31 

2 

0.438, 0.582 

Nonlinear Optimization 

0.14 

12.98 

-0.71 

0.35 

3 

0.431, 0.917, 1.314 

Nonlinear Optimization 

0.07 

19.22 

0.58 

-0.01 

4 

0.511, 0.802, 0.820, 1.505 

Nonlinear Optimization 

0.06 

19.87 


0.09 


In order to determine sets of aerodynamic 
approximations which consistently improve a chosen 
criterion as the order of the approximation increases (with 
increasing number of denominator coefficients), nonlinear 
optimization of the denominator coefficients must be 
employed. In this case the criterion for goodness of fit is 
the total physically weighted error. Table 4 shows the 
comparison of total physically weighted error ( w cut = 0), 
as well as the total normalized error of refs. 9-11, and % 
error in flutter q and frequency for four optimized, 
physically weighted MS approximations. Figure 5 shows 
the decreasing total physically weighted error criterion 
with increasing order. As would be expected since the 
selection of denominator coefficients is based on optimum 
physically weighted errors, the normalized error is poor 
and does not improve with increasing order. There is no 
corresponding improvement in flutter q and frequency, but 
figure 6 shows that there is an improving trend in overall 


frequency and damping characteristics at flutter of all the 
elastic modes with increasing order as compared to the 
frequency domain baseline flutter root characteristics. 



Figure 5.- Comparison of total weighted errors for opti- 
mized physically weighted MS approximations 
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Figure 6.- Comparison of open loop flutter root characteristics for optimized physically weighted approximations 


Imag 

(rad/sec) 



(0 < q < 400 psf) 

%difference in flutter q < .8 (=2.6 psf) 
%difference in flutter frequency < .4 (=.3 rad/sec) 


Figure 7.- Stability analysis root loci for the AFW antisymmetric modes at Mach 0.9 using optimized physically 
weighted Minimum-Slate approximations. 


Figure 7 shows the stability analysis root loci of 
the 1-state and 2-state optimized, physically weighted MS 
approximations. The figure demonstrates that the 
optimized 1 -state case matches the baseline case fairly 
well and that extremely close correlation is achieved with 
the optimized 2-state approximation. The statements 
made with respect to the roll mode root locus in 
discussion of figure 4 apply here as well. 


Figures 6 and 7 demonstrate that the physical 
weighting criterion for the vibrational modes is effective 
in achieving approximations which yield accurate stability 
characteristics for poles with low damping-ratios. It is 
interesting that agreement with the p-k root loci are also 
good for relatively large damping ratios despite the fact 
that the p-k approximation Q(p) = Q(0+ik) is degrading; 
one only finds in the case of the roll mode where the p-k 
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approximation is clearly invalid that there is no consistent 
improvement in damping ratio as compared to the baseline 
with increasing order of approximation. Results (not 
shown) of Minimum-State approximations with the same 
equality constraints shown herein which reduce normalized 
error have higher % errors in flutter q and frequency and are 
much more sensitive to optimization of lag coefficients. 
Physical weighting tends to improve open-loop flutter 
characteristics and reduces the sensitvity to choices of 
aerodynamic roots. Corresponding control and gust 
response analyses, as well as closed loop flutter analyses 
would have to be presented in order to demonstrate the 
goodness of the physical weighting criteria for the control 
and gust modes. 


By combining various capabilities, namely the 
Minimum-State method, selectable constraints, optimal 
selection of denominator coefficients, as well as the 
determination and use of critical frequency ranges for 
approximating the tabular generalized aerodynamic forces 
with rational functions based upon physical properties, it 
is possible to obtain good, low order state-space 
mathematical models for design and simulation of 
aeroservoelastic systems. A state-space model with only 
two aerodynamic states predicted a flutter dynamic pressure 
with less than 0.8 percent error and showed good 
agreement in the root loci of all elastic modes as compared 
to conventional stability analysis. The significance of 
these results is that good models can be developed having 
an order of magnitude fewer augmenting aerodynamic 
equations than more traditional approaches. This reduction 
facilitates the design of lower order control systems, 
analysis of control system performance, and near real-time 
simulation of aeroservoelastic phenomena. 
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